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MISS  DISTANCE  MEASUREMENT  USING  ONE  CAMERA 

by 

A.  T.  Smith 

SUMMARY 

Data  from  one  high  frame  speed  camera  on  a  target  can  give  missile 
trajectory  and  miss  distance.  This  Report  gives  a  theoretical  solution  and 
indicates  how  it  can  be  modified  to  process  practical  data  in  a  computer.  The 
sources  of  error  and  their  magnitudes  are  given. 
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INTRODUCTION 


Whenever  a  missile  is  fired  at  a  target  it  is  important  during  the  final 
stage  of  attack,  to  obtain  relative  trajectory  and  miss  distance  information. 

At  Aberporth  these  are  either  aerial  or  sea  targets  and  the  attacking  missile  is 
normally  sighted  by  one  or  more  high  speed  cameras.  These  cameras  are  free 
running  and  run  nominally  at  one  hundred  frames  per  second.  Time  information 
from  the  Range  Central  Timing  Unit  is  recorded  on  the  edge  of  the  camera  film, 
so  that  an  accurate  time  can  be  given  to  each  frame. 

Only  two  types  of  aerial  targets  at  present  carry  a  camera  system,  these 
are  aircraft  targets  and  towed  targets  (infra-red  or  radar  sources  towed  by  an 
aircraft).  Normally  the  missile  is  photographed  by  two  camera  packs  (each  pack 
consisting  of  two  cameras)  one  on  each  wing  tip  of  the  target  aircraft.  The  data 
obtained  from  these  cameras  can  then  be  interpolated  to  corresponding  times,  to 
allow  a  simple  two  camera  solution^  using  the  wing  span  of  the  aircraft  as  a 
baseline.  It  has  been  found  from  experience  that  the  extrapolation  of  camera 
data  seldom  gives  satisfactory  results. 

For  aircraft  targets  and  13  metres  towed  targets  the  camera  packs 
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consist  of  two  WRETAR  Mark  3  cameras  .  These  cameras  have  a  field  of  view  of 
186°  and  are  mounted  to  give  a  complete  spherical  optical  coverage.  At  least 
three  occasions  arise  on  these  targets  when  a  single  camera  solution  must  give 
trajectory  d'ata.  These  occur  when 

(a)  a  faulty  camera  fails  to  record  either  pictures  or  the  timing  information 

(b)  a  camera  pack  is  removed  to  allow  additional  equipment  on  the  aircraft, 

(e.g.  a  telemetry  pack) 

(c)  the  time  interval  common  to  both  cameras  for  interpolating  data  is  too 

short  to  allow  a  satisfactory  two  camera  solution  over  the  intercept  period. 

Although  the  camera  packs  have  full  spherical  coverage,  on  some  flight  paths  the 
missile  is  obscured  by  the  aircraft,  for  part  of  the  intercept,  from  one  or 
both  cameras.  This  situation  can  be  made  worse  by  certain  events  (e.g.  flash) 
which  results  in  the  loss  of  further  data.  Quite  often  in  these  conditions  one 
or  both  cameras  will  give  sufficient  data  for  a  single  camera  solution,  given 

an  independent  velocity  estimate. 

The  single  camera  solution  is  considered  too  inaccurate  for  use  on  the 
longer  range,  30  metres  and  61  metres,  towed  targets.  The  greater  ranges 
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involved,  of  the  missile  from  the  camera,  (see  section  5)  and  the  probable 
inaccuracies  in  the  knowledge  of  the  tow  position  are  the  main  reasons  for  this 
decision. 

The  Range's  sea  targets  were  basically  designed  for  a  single  camera 
solution  as  the  fields  of  view  of  the  target  cameras  do  not  overlap  to  any 
great  extent,  except  on  the  largest  targets.  The  cameras  used  are  GW1  Mark  1A 
and  GW2  Mark  2,  and  these  have  a  field  of  view  of  120°.  This  means  that  most 
missile  trials  against  sea  targets  require  a  single  camera  solution. 

Until  recently  a  manual  single  camera  solution  described  in  a  previous 
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Report  has  been  used.  This  manual  solution  has  now  been  superseded  by  a 
computer  solution.  This  Report  describes  the  theoretical  solution  of  the 
problem  and  how  this  solution  is  modified  for  computer  use  in  the  practical 
case. 

2  THEORETICAL  SOLUTION 

To  obtain  a  theoretical  solution  for  trajectory  data  using  a  single 
camera,  certain  assumptions  have  to  be  made  about  the  trajectory  during  the 
period  of  time  covered  by  the  solution.  The  trajectory  in  question  is  the 
relative  trajectory  of  the  missile  with  respect  to  (a)  the  camera  origin  and 
(b)  the  camera  axes  system.  These  assumptions  are 

(1)  the  missile  flies  in  a  straight  line 

(2)  the  velocity  V  of  the  missile  is  a  known  constant. 

In  practice  these  assumptions  imply  that  the  target  (and  hence  camera)  is 
also  moving  in  a  straight  line  at  constant  velocity  and  that  the  attitude  of 
the  target  is  steady. 

Since  basically  any  camera  system  supplies  direction  cosine  vectors  of  a 
trajectory  with  respect  to  the  camera,  the  following  theory  is  based  on  direc¬ 
tion  cosine  vectors.  For  the  theoretical  solution,  it  will  be  assumed  that 
there  are  no  errors  in  the  direction  cosine  vectors.  The  effect  of  these  and 
other  practical  errors  will  be  discussed  briefly  in  a  later  section. 

The  missile  flies  in  a  straight  line  and  therefore,  provided  the  camera 
does  not  lie  in  the  path  of  the  missile  the  missile  trajectory  and  the  camera 
origin  define  a  unique  plane  in  space.  In  addition  all  rays  from  the  camera  to 
the  missile  and  hence  all  direction  cosine  vectors  of  the  missile  with  respect 
to  the  camera,  lie  in  this  plane.  The  missile  also  flies  at  a  constant  velocity 
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both  spatially  and  with  respect  to  camera  axes,  and  therefore  the  distance  flown 
between  two  given  time  points  depends  only  on  the  time  difference  and  not  on 
absolute  time. 

Let  P.^  =  (x£,  y^,  z^)  be  the  position  vector  of  the  missile  with  respect 
to  the  camera,  at  the  time  t^  .  Then  a  constant  velocity  implies  that 

D.  <=  t.  ,  -  t.  , 

1  l+l  l  * 


t.  ,  -  t. 
l+l  i 


=  constant  (  =  velocity  of  the  missile) 


where  D.  =  P 


■i+1  -i' 


■JK+ 1  -  xi>2  +  <yi*i  -  yi>2  *  (zi+i  -  2P2} 


see  Fig. 1 . 


2 . 1  Calculation  of  missile  direction 

Consider  any  line  parallel  to  the  missile  trajectory  which  also  lies  in 
the  camera-trajectory  plane.  The  length  of  a  section  of  this  line,  formed  by 
two  rays  from  the  camera  to  the  missile  trajectory,  depends  only  on  the  time 
difference  between  these  two  rays.  This  follows  because  the  missile  trajectory 
also  has  this  property.  Therefore  the  direction  of  the  missile  can  be  found 
from  just  direction  cosine  vectors.  This  is  achieved  by  finding  a  line  in  the 
camera- trajectory  plane  with  the  following  property :- 

The  direction  cosine  vectors  from  the  camera  should  cut  the  line  into 
sections  whose  length  depend  only  on  the  time  difference  between  the  successive 
direction  cosine  vectors. 

With  reference  to  Fig. 2  let  cL  be  the  direction  cosine  vector  of  the 
missile  at  time  t^  .  Let  be  the  length  of  the  section  of  a  line  formed  by 


the  successive  direction  cosine  vectors  d.  , 

—x+] 


and  d.  .  Then 

—l 


fci+l 


=  constant  for  a  line  parallel  to  the  missile  trajectory 


A  line  with  this  property  can  be  obtained  by  using  the  following  procedure 
Define  a  new  system  of  axes  X’,  Y’ ,  Z'  such  that 
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(a)  Z'  is  in  the  direction  of  a  vector  perpendicular  to  the  camera-trajectory 
plane , 

(b)  Y'  is  in  the  direction  of  the  first  direction  cosine  vector  dj  , 

(c)  X'  is  in  the  direction  of  a  vector  perpendicular  to  both  the  Y*  and  Z* 
axes,  and 

(d)  X',  Y'  and  Z'  form  a  right  handed  triad. 

Both  the  new  X'  and  Y'  axes  lie  in  the  camera-trajectory  plane  and 
hence  form  a  base  for  the  plane.  Any  direction  cosine  vector  in  this  plane  can 
now  be  described  by  the  single  variable 


y  .  d. 

S.  =  - -L  (see  Tig. 3) 

where  x  and  y  are  unit  vectors  along  the  axes  X?  and  Y*  respectively. 
Note  that  y  =  S^x  is  a  line  in  the  plane  (X',  Y')  which  lies  in  the  same 
direction  as  d_^  and  that  S  j  =  °°  . 

Since  any  line  with  the  relevant  property  is  required,  a  line  of  the  form 

y  =  Mx  +  1 

is  calculated.  This  line  intersects 


when 


y  =  S.x 
i 


S.x  =  Mx  +  1  , 

l  ’ 


i.e. 


x 


and 


y 


In  theory  three  direction  cosine  vectors  are  now  required.  For  convenience  let 
these  be  the  first  three  vectors.  The  line  y  =  Mx  +  1  intersects  the  lines 
y  =  SjX,  y  =  S£X  and  y  =  S^x  at  the  points 
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(Xj,  y,)  =  (0,  1) 


^x2 *  y2^  fs2  -  M  *  S2  -  M  +  1 


(x3’  y3*  (s,  -  M  *  S0  -  M 


+  3  I  respectively. 


See  Fig. 4.  Then 


\(  1  l  V  .  [  m  mV 
2  n  1\s3  -  M  S2  -  M y  l  S3  -  m  s2  _  My 


L  .  m2n  _ 1 _ i_ 

(1  M  ^  S3  -  M  S2  -  M 


For  the  line  y  =  Mx  +  1  to  have  the  required  property 


t2  tj  t3  -  t2 


where  tj,  t2  and  t3  are  the  times  of  the  first  three  direction  cosine 


vectors. 


Therefore 


(S2  -  M)(t 


(t2  -  tj)  c3  ~  t2  (^3  ”  M  S2  ”  M} 


therefore 


1  r  3  "  fcl 


S2  "  M  \t2  ”  fcl 


S3  -  M 


therefore 
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Hence 


s3(t3  -  t,)  -  s2(t2  -  tj) 


t3  t2 


Note  that  the  x  and  y  co-ordinates  of  the  missile  in  the  (X',  Y')  plane 
are  both  linear  in  time  t  .  It  is  therefore  possible  to  find  a  pair  of 
equations  of  the  form 

y  =  m(t  -  tj)  +  1 


x  =  n (t  -  tj)  , 


instead  of  the  single  equation  y  =  Mx  +  1  .  The  auxiliary  equation  y 
then  automatically  has  the  required  property. 


Therefore 


Q  ^ 

S2  n 


S  -  - 
3  n 


Therefore 


S2  S3 


1  1 


n  2  "  fcl  t3  *1 


therefore 


_ fc3  t2 _ 

(s2  -  s3)(t2  -  tj)(t3  -  tj) 


S3(t3  -  tj)  -  S2(t2  -  tj) 
m  =  (S2  -  S3)(t2  -  tj)(t3  -  tj)  • 

In  the  computer  program  solution  of  the  problem,  the  solution  is  simplified  by 
calculating  m  and  n  . 

2.2  Calculation  of  missile  trajectory 

Now  the  trajectory  of  the  missile  has  an  equation  of  the  form 


y  =  n  x+  r  * 


3  13 


in  the  (X1,  Y’)  plane  and  the  co-ordinates  have  equations  of  the  form 


y  =  m(t  -  tj)r  +  r 
x  =  n(t  -  tj)r 


- 2 - [n(t  “  t])£+  {m(t  “  tj)  +  1  }y] 

/(m  +  n  ) 


This  now  completes  the  theoretical  solution  since  it  is  possible  to  calculate  the 
missile  position  at  any  given  time.  The  calculation  of  the  time  and  range  of 
nearest  miss,  and  other  useful  parameters  is  given  in  section  4. 

3  PRACTICAL  SOLUTION 

Considering  a  missile  travelling  at  700  m/s,  with  a  limit  of  25  g  on  the 
lateral  acceleration,  it  can  be  shown  that  there  exists  a  straight  line  from 
which  the  missile  cannot  deviate  by  more  than  16  cm,  over  a  period  of  1/10  s 
(a  likely  maximum  time  for  a  missile  to  be  in  view).  Therefore  in  practice,  a 
missile  will  not  deviate  much  from  the  assumed  model,  at  least  while  passing 
the  camera.  This  model  can  therefore  be  used  as  an  approximation,  to  obtain  a 
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trajectory.  The  solution  requires  the  velocities  of  the  missile  and  the  target, 
which  must  be  supplied  from  external  sources.  Due  to  errors  introduced  by  the 
camera  system  and  the  error  of  applying  the  model  to  the  missile,  the  solution 
needs  to  be  modified  in  the  practical  case. 

The  following  information  will  be  assumed  for  the  solution  of  a  practical 
situation : 

(a)  At  least  three  direction  cosine  vectors  of  the  missile  with  respect 
to  a  camera  and  their  corresponding  times,  i.e. 

(ti ,  d i)  for  i  =  1 ,2 , . .  . ,r (r  >  3) 

where  d.  =  (£. ,  m. ,  n.)  ,  and 
— l  i  l  l 

(b)  the  relative  velocity  V  of  the  missile  or  at  least  an  estimate  for 
it. 

3. 1  Calculation  of  the  camera-trajectory  plane 

The  first  stage  of  the  computer  solution  is  the  calculation  of  the  camera- 
trajectory  plane.  Since  the  direction  cosine  vectors  will  not  lie  perfectly 
in  a  common  plane,  the  'best'  plane  passing  through  them  must  be  defined  and 
then  calculated.  Before  continuing  it  must  be  noted  that  any  plane  passing 
through  the  origin  can  be  defined  by  a  single  unit  vector  (L,  M,  N)  .  This 
vector  is  perpendicular  to  any  vector  lying  in  that  plane.  If  (L,  M,  N) 
satisfies  this  condition  then  so  does  (-L,  -M,  -N)  ,  therefore  it  will  be 
assumed  that  N  >  0  .  With  reference  to  Fig. 5  define  0^  by 

LI.  +  Mm.  +  Nn.  =  cos  [5  -  0.  } 

1  1  1  V2  V 

=  sin  0. 

l 

The  'best'  plane  is  defined  as  that  plane  which  minimizes 


r 


i=  1 


=  ^  '  |sin  '  +  Mnu  +  NnJ  | 

i=  1 


S(L,M,N) 
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Now  for  small  0  ,  sin  0^0^  sin  *  0  .  Therefore  the  problem  is  simplified  by 
finding  the  vector  (L,  M,  N)  which  minimizes 


S(L,M,N) 


*  2 
X  {L£i  +  ^i  +  Nni} 

i=l 


2  2  2 

Now  since  L  +  M  +  N  =1  and  N  is  assumed  non-negative  the  problem  is  to 

find,  say  L  and  M  ,  such  that  S(L,M)  is  a  minimum,  where 


S  (L,M) 


+  Mm. 

l 


S(L,  M)  has  a  minimum  when 


BS 

3L 


3S. 

3M 


0 


y 


i= 1  i=l 


12  4  - 12  1  *  m  12  * --  v k  *  E  ni’-i  -  f-  12  mi°t  ■ 


i=  1 


i=  1 


i=  1 


and 


M 


i=  1  i=  1  i=l  i=l  i=  1 
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The  root  of  the  above  two  non-linear  equations  gives  the  required  L  and  M  . 

The  problem  of  finding  the  best  plane  has  now  reduced  to  solving  these  two  non- 

.  2  2  2  2 

linear  equations.  Since  -  £n^,  Ef-nu,  £n^£^,  Enun^  and  Em.  -  En.  are 

constants  they  will  be  replaced  by  a,  b,  c,  d  and  e  respectively. 

Let 


f , (M,L) 


La  +  Mb  + 


(N2  -  L2)c 
N 


ML  , 
d 

~T 


o 


and 


f 2 (M,L) 


Me  +  Lb  + 


(N2  -  M2)d 
N 


0  . 


The  root  of  this  pair  of  non-linear  equations  is  found  by  using  the  Newton- 
Raphson  technique.  The  details  of  this  technique  are  given  in  Appendix  A.  Let 
(M,  L)  be  the  root  required  and  (M,  L)  be  an  approximation  to  this  root,  then 


and 


M  *  M  - 

gy  -  a 


L  *  L  -  ~  aA 

$Y  ~  a 


where  a  = 


6  = 


3f  j  (M,L)  3f2(M,L) 


Y  = 


3M 

3f, (M,L) 
3L 

3f 2 (M,L) 


3L 


3M 

X  =  f  (M,L)  and 
4)  =  f2(M,L)  . 


3f, (M,L) 
3M 

3  f  j (M,L) 
3L 


b  +  ~  |(M2  -  1 )Mc  +  (L2  -  l)Ld| 


=  a  + 


N' 


|(3M2  +  2L2  -  3)Lc  +  (M2  -  l)Md| 
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3f2(M,L) 

3M 


e  + 


n'3  ^ 


(3L2  +  2M2  -  3)Md  +  (L2 


3f2(M,L) 

3L 


b  + 


4{< 

NJ  1 


(M  -  1 )Mc  +  (L  -  1 )Ld 


}• 


The  Newton-Raphson  technique  gives  rise  to  an  iterative  process  for  finding 
the  root  (L,  M)  .  An  initial  estimate  for  the  root  can  be  obtained  by 
considering  the  plane  formed  by  two  of  the  original  direction  cosine  vectors. 
Using  the  first  and  last  direction  cosine  vectors,  the  initial  estimate  for  L 
and  M  is  given  by 


(L,  M,  N) 


if  £,m  -  m,£  ^  0 

1  r  1  r 


if  £.m  -  m.£  <  0 

1  r  1  r 


where  d.  =  (£.  ,  m. ,  n. ) 

—l  l  l  l 

and  d.  *  d  =  (m,n  -  n,m  ,  n.5,  -  n  £,m  -  m.H  ) 

—1  — r  1  r  1  r  1  r  r  1 ’  1  r  1  r 

Note  both  functions  fj(M,L)  and  f2(M,L)  and  their  derivatives  have 
poles  at  N  =  0  ,  (i.e.  =  1 )  .  Therefore  if  N  is  small  or  zero,  the 

iteration  may  not  converge.  These  poles  have  no  physical  significance  and  are 
a  consequence  of  the  axes  system  used.  Therefore  if  N  is  small  or  zero  a 
change  to  another  system  of  axes  X' ,  Y',  Z'  ,  where  N'  is  not  small,  allows 
the  vector  (L',  M' ,  N')  to  be  calculated.  This  by-passes  the  poles  and  the 
problem  of  convergence.  Once  the  vector  (L* ,  M' ,  N')  has  been  calculated  a 
rotation  back  to  the  original  system  of  axes  gives  (L,  M,  N)  . 


At  this  stage  the  'best'  plane  nearly  containing  the  direction  cosine 
vectors  ck  has  been  calculated.  Now  in  general  these  direction  cosine  vectors 
will  not  lie  in  this  plane,  therefore  these  vectors  must  be  projected  onto  the 
plane.  Let  £  =  (L,  M,  N)  ,  then  the  projection  of  cL  onto  the  camera- 
trajectory  plane  is 


-  (P  .  dj.JP  . 
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This  vector  will  not  in  general  be  a  unit  vector,  therefore  it  should  be  normal¬ 
ized  before  continuing.  Let  e^  be  this  modified  direction  cosine  vector,  then 

d.  -  (P  .  d.)P 

—l  —  — i  — 

e .  =  - 

Id.  -  (P  .  d.)p| 

3. 2  Calculation  of  missile  direction 

Define  a  new  system  of  axes  X',  Y'  and  Z' 
where  X*  is  in  the  direction  of  e.  *  P_  =  Q  say 

Y'  is  in  the  direction  of 
and  Z'  is  in  the  direction  of  P^  =  (L,  M,  N)  . 

Q  =  (g,N  "  hjM,  h jL  -  Nfp  f,M  -  g,L)  , 
where  e.,  =  (f£ ,  g£,  h. )  . 

This  gives  a  right  handed  system  of  axes  Q,  e^  and  P  .  Using  the 
vectors  Q  and  which  both  lie  in  the  camera-trajectory  plane,  define  the 

variable  S  at  the  time  t 


S. 

i 


Note  Sj  =  °° 

At  this  stage,  the  best  pair  of  equations  of  the  form 


y  =  m(t  -  tj)  +  1  and  x  =  n(t  -  tj) 

to  satisfy  the  assumed  model,  have  to  be  defined  and  then  calculated.  The  curve 
y  =  mx/n  +  1  intersects  the  lines  y  =  S.x  when 


nS. 

l 


m 


and 


y  = 


m 


nS.  -  m 

l 


+  1 


In  theory  these  co-ordinates  should  be 


x 


n(ti  ■  V 


and 


y 


m(t^  -  t  j )  +  1 
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Therefore  the  best  values  of  m  and  n 


are  defined  as  those  which  minimize 


A(m 


•“>  '  12  '  V  '  Sspu)  +  (n(ti  '  V  "  5sfr 


mi 


(m2  + 


-2>E(h-h-^)2  • 

i=2 


A(m,n)  has  a  minimum  when  =  0  ,  i.e. 

v  *  '  3n  Sm  ’ 


12  (£i  •  h  -  '  (m2  +  “2)  12  (h  '  h  •  5s p-s) 


m  y  It.  -  t 
i=2 

and 


1 


1  nS-  •  m'  <„s.  -  „)2 


=  0 


n  V  ft.  -  t,  -  <_f  .  (m2  *  n2)  V  ft.  -  t,  -  -^±-)  — !i  , 

££  V  1  1  nSi  -  ”7  ^  \  1  1  nS1  -  "/  (ns.  -  m>2 


=  0 


The  root  of  this  pair  of  non-linear  equations  gives  the  best  values  of  m  and 
n  .  Again  the  Newton-Raphson  technique  is  used  to  find  the  root. 

Let 

gj (m,n) 


and 

g2 (m,n) 


Let  (m,  n)  be  the  root  and  (m,  n)  be  an  approximation  to  the  root,  then 


m 


m 


a<j>  -  yX 
2 

a  -  By 


and 


n  ^  n 


otX  -  $<j> 

2  * 
a  -  By 
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where  a  = 


6  = 


3g  (&,ii)  3g  (m,n) 


Y  = 


3m 

3gj(m,n) 

3m 

3g2(m,n) 

9n 


3n 


X  =  g  j  (m,fi)  and 
<j>  =  g2(m,h) 


3gj (m,n) 
3n 


+  2(m2  +  n2)  7;  (t.  -  t, 
i=2  ' 


S. 

_ l 

(nS^  -  m) 


3 


(m2  +  n2)  ^ 
i-2 


3g2 (m,n) 
3m 


3  g  j (m»n) 
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An  initial  estimate  for  this  iteration  can  be  obtained  using  the  data  of  three 
time  points.  For  example,  tj,  t2  and  t^,  then 


and 


<S2  S3^  ^2  ~  tl^t3 


s3(t3  t})  -  s2(t2  -  tj) 

(S2  -  S3)(t2  -  tj)(t3  -  tj) 


Calculation  of  missile  trajectory 


The  missile  trajectory  is  now  calculated  as  in  the  theoretical  solution. 
The  position  vector  of  the  missile  at  the  time  t  ,  in  camera  axes  and  with 
respect  to  the  camera  origin  is  therefore 


-  v  +  • 


The  target  origin  and  axes  system  is  in  general  different  to  that  of  the  camera, 
therefore  to  complete  the  solution  the  trajectory  has  to  be  converted  to  this 
new  origin  and  axes  system.  The  actual  target  origin  is  normally  a  nominal 
position  on  the  target,  for  example  the  centre  of  gravity  of  the  aircraft. 

This  solution  has  been  written  into  a  computer  program  subroutine,  a 
listing  and  short  description  is  given  in  Appendix  B. 

4  HISS  DISTANCE  PARAMETERS 

This  section  gives  the  mathematics  for  the  calculation  of  the  normal 
parameters  required  at  the  time  of  nearest  miss.  The  time  of  nearest  miss  is 
the  time  at  which  the  range  of  the  missile  from  the  target  is  a  minimum.  The 
parameters  required  at  nearest  miss  are 
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(a)  Time  t 


miss 


(b)  Range  R  . 

°  miss 

(c)  Co-ordinates  (xmiss>  Ymiss’  Zmiss^ 

(d)  Direction  cosine  vector  of  the  trajectory  (&,  m,  n)  and 

(e)  Velocity,  assumed  supplied  from  other  sources. 

Let  the  position  of  the  missile,  with  respect  to  the  target,  at  time  t  , 
be  (xj,  y j ,  Zj)  and  at  time  t2  be  (x2>  y2,  z£)  .  Then  the  position  of 
the  missile  at  time  t  is 

(xn  +  Lt,  yn  +  Mt,  z  +  Nt) 


1 


where  xQ  =  x,  +  ^  -  t] 

fcl 


yo  =  yl  +  t  -  t 
u  c2  1 

C1 


z0  *  Z1  *  4  -  t, 


t2x,  -  tfy 
<XI  -  x2>  -  t-  -  t'  — 


(y,  -  y,)  - 


1 


t2y  1  ~  fcly2 
t2  "  t] 


s  t2Z 1  tlZ2 

(zl  "  Z2  "  t2  -  tj 


L  = 


x2  ~  X1 

t  —  t  ’ 

2  1 


M  = 


y2  ~  yi 
fc2  ”  C1 


and  N 


z2  -  Z1 
t2  “  fcl 


L,  M  and  N  are  velocity  components. 
The  velocity  of  the  missile  is 


V 


J{L* +  «2  +  n2}  ■  irrh ;J{(X2  ~  xi)2  *  <y2  -  V2  +  <z2  -  zi>2} 


and  the  direction  cosine  vector  of  the  missile  trajectory  is  given  by 


(£,  m,  n)  =  -(L,M,N) 


1 


VI 


(x2  -  Xj)2  +  (y2  -  y,)2  +  (z2  -  z^2 


(x2  -  ,  y2  -  yj,  z2  "  ZP  * 


Now  the  range  of  the  missile  from  the  target  is  given  by 
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There  are  three  sources  of  error  for  the  single  camera  solution,  in  a 
practical  situation,  these  are:- 

(a)  Reading  errors 

Both  time  and  angular  data,  have  to  be  read  from  films.  This  can  only 
be  done  to  a  finite  accuracy  which  depends  on  the  camera  system  used. 

The  timing  information  on  the  Range's  films  can  be  read  to  an  accuracy  of 
0. 1  milliseconds,  therefore  the  maximum  percentage  error  in  the  calculated 
trajectory  data,  due  to  a  timing  error  is 


_  0.0001 

2  x  — -  x 


100% 


0  ■  02o, 

T 


(T  in  seconds) 


Let  a  be  the  accuracy  of  the  angular  data  read  from  the  film  in  radians, 
then  the  maximum  percentage  error  is  100a%  which  is  probably  pessimistic  since 
the  least  squares  technique  used  should  reduce  this  error. 

(b)  The  error  in  applying  the  model,  required  by  the  single  camera 
solution,  to  an  actual  trajectory 

The  missile  can  fail  to  satisfy  the  model  for  at  least  three  reasons 

2 

(i)  a  lateral  acceleration,  of  say  'a'  m/s 

2 

(ii)  a  longitudinal  acceleration,  of  say  m/s  ,  and 

(iii)  a  rotating  camera  axes  system,  say  the  axes  system  is 
rotating  at  a  rate  of  'b'  rad/s. 

(N.B.  The  lateral  and  longitudinal  accelerations  are  in  general  the  combined 
effects  of  both  the  missile  and  the  camera  moving.) 

It  can  easily  be  shown  that  the  maximum  deviations  from  the  model  are 


(D 

(2) 

and 


metres,  due  to  the  lateral  acceleration 


metres,  due  to  the  longitudinal  acceleration 


(3) 


2 

bVT 

8 


metres,  due  to  the  rotating  axes  system. 
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(c)  The  error  in  the  velocity  estimate 

A  5%  error  in  the  velocity  will  produce  a  5%  error  in  the  calculated 
trajectory  data.  (It  should  be  noted,  that  all  the  percentage  errors  quoted 
above,  apply  only  when  the  trajectory  origin  is  identical  to  the  camera  origin.) 

Given  the  characteristics  (actual  or  theoretical)  of  the  missile,  target 
and  camera  system  (e.g.  velocity,  lateral  accelerations,  etc.),  it  is  therefore 
possible  to  calculate  the  maximum  error  to  be  expected  from  using  a  single 
camera  solution.  For  example,  consider  the  representative  figures  for  a  missile 
fired  against  an  aerial  target,  typical  maximum  values  would  be 

R  =  35  m 

a  =  150  m/s2  («  15  g) 

£  =  50  m/s2  («*  5  g) 

b  =  10°/s 

i  o 

a  =  £ 

V  =  700  m/s  with  a  5%  error. 

Then  the  maximum  compounded  error  in  position,  assuming  a  typical  time  T 
of  0.05  s,  is  approximately  2.27  m,  of  which  1.75  m  is  due  to  the  velocity  error. 
It  should  be  noted  that  this  compounded  error  has  its  maximum  value  at  the 
extreme  range  from  the  camera  and  will  be  considerably  less  nearer  the  camera. 

The  main  source  of  error  for  the  single  camera  solution,  in  the  example 
quoted,  is  the  velocity  estimate.  This  is  considered  to  be  the  situation  for 
most  single  camera  solutions  calculated  at  Aberporth.  The  reason  being  that  the 
tracking  radars  usually  fail  to  give  trajectory  data  and  hence  velocity  data 
during  the  period  of  the  intercept. 

When  the  trajectory  data  is  calculated  with  respect  to  an  origin  and 
axes  system,  other  than  that  of  the  camera  (e.g.  if  the  target  aircraft  is 
towing  a  decoy  target  behind  it),  additional  errors  can  easily  be  introduced. 
These  errors  are  obviously  dependent  on  the  accuracy  of  the  trajectory  origin 
and  axes  system  with  respect  to  those  of  the  camera,  but  must  be  smaller  for 
origins  nearer  the  camera. 

6  CONCLUSION 

The  single  camera  solution  puts  fairly  restrictive  conditions  on  the 
trajectory  of  the  missile.  These  conditions  limit  the  situations  in  which  a 
single  camera  solution  can  be  used.  This  is  in  contrast  with  solutions 
obtained  from  two  or  more  cameras.  For  example,  the  single  camera  solution 


22 


162 


is  unsuitable  for  obtaining  trajectory  data  during  the  launch  phase  of  a  missile. 
This  is  because  the  missile  is  undergoing  substantial  longitudinal  acceleration 
at  relatively  low  speed,  during  this  period. 

Another  restriction  on  the  single  camera  solution  is  the  fact  that  the 
velocity  of  the  missile  relative  to  the  camera  must  be  supplied  from  an  external 
source. 

Nevertheless,  taking  these  restrictions  into  account  the  single  camera 
solution  can  be  very  useful  in  obtaining  trajectory  data.  This  is  particularly 
so  during  the  attack  phase  of  a  missile,  when  the  missile  trajectory  normally 
approximates  to  the  model  required  by  the  single  camera  solution.  The  solution 
is  not  specifically  restricted  to  aerial  or  sea  targets  and  will  apply  to  any 
single  camera  situation  where  the  necessary  extra  information  (timing,  missile 
velocity,  etc.)  is  available. 
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Appendix  A 

DETAILS  OF  THE  * NEWTON-RAPHSON *  TECHNIQUE 

Given  here  are  the  details  for  finding  the  root  of  a  pair  of  non-linear 
equations,  using  the  Newton-Raphson  technique. 

Let  the  pair  of  non-linear  equations  be 


and 


fj(M,L)  =  0 


f2(M,L)  =  0  . 


Let  (M,  L)  be  the  root  and  (M,  L)  be  an  approximation  to  this  root,  then 


f .  (M,L)  «  f.  (M,L)  -  6M 


3f.(M,L)  3f . (M,L) 


3M 


-  6L 


3L 


and 


3f 2 (M,L)  3f2(M,L) 


f2(M,L)  *  f2(M,L)  -  6M  -  -  6L 


9L 


where  6M  =  M  -  M  and  6L  =  L  -  L 
Since 

f  j  (M,L)  =  f,(M,L)  =  0 

it  follows  that 


and 


Let 


9f  (M,L)  3f  (M,L) 

*  •  Pt  * 


f,(M,L)  »  <SM  - - -  +  6L 


3L 


3f  (M,L)  3f „ (M,L) 

f2(M,L)  - 


3fj(M,L) 

3M 

3f2(M,L) 

3L 


=  a 


=  to 


9fj(M,L) 


3L 


=  6 


3f2 (M,L) 


3M 


f  j  (M,L)  =  X  and  f2(M,L)  =  *  , 
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then 


A  »  a6M  +  BfiL 


and 


<j)  ^  y6M  +  w6L 


Therefore 

SM  «  |.t- 2X 
By  ~  wa 


and 


sl  „  li-iai 

By  -  wa 


therefore 


M  «  M  -  6M 


and 


L  L  —  6L 


This  gives  rise  to  an  iterative  process  for  calculating  the  root  (M,  L)  . 
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Appendix  B 

SHORT  DESCRIPTION  AND  LISTING  OF  THE  SINGLE  CAMERA  SOLUTION  SUBROUTINE  SINCA 

The  subroutine  SINCA  performs  the  single  camera  solution  described  in 
Section  3.  Starting  with  a  velocity  estimate  and  at  least  three  direction 
cosine  vectors  and  their  corresponding  times,  the  subroutine  calculates 
trajectory  data  and  miss  distance  information  with  respect  to  any  given  origin 
and  axes  system.  Trajectory  data  is  automatically  produced  at  the  times  of  the 
direction  cosine  vectors  and  the  time  of  nearest  miss.  Since  trajectory  data  is 
often  required  at  special  times,  e.g.  at  fuse  triggering  time,  there  is  an 
option  to  supply  data  at  other  given  times.  The  subroutine  is  called 

SINCA  (T,  RL,  RM,  RN,  IN,  JN,  X,  V,  RMISS,  D) 

T,  RL,  TM,  RN,  X,  RMISS  and  D  are  arrays  dimensioned  as  follows:  T(25),  RL(25), 
RM(25) ,  RN (25) ,  X(3),  RMISS(8)  and  D(3,3). 

INPUT  DATA 

IN  :  number  of  direction  cosine  vectors  (3  <  IN  <  24) 

JN  :  number  of  special  time  points  (IN+JN  *5  25) 

T(I)  :  (a)  time  of  Ith  direction  cosine  vector  I  =  1 ,  IN 
(b)  time  of  a  special  time  point  I  =  IN+1,  IN+JN 

RL(I)  :  1  component  of  Ith  direction  cosine  vector  I  =  1,  IN 

RM(I)  :  m  component  of  Ith  direction  cosine  vector  I  =  1 ,  IN 

RN(I)  :  n  component  of  Ith  direction  cosine  vector  I  =  1,  IN 

X  :  (x,  y,  z)  co-ordinates  of  the  target,  with  respect  to  the  axes  system  of  the 
direction  cosine  vectors  and  the  camera  origin. 

D  :  rotation  matrix 

V  :  velocity  estimate 
OUTPUT  DATA 

IN  :  as  input 
JN  :  as  input 

V  :  as  input 
X  :  as  input 
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D  :  as  input 


T  (I) 

:  as 

input  I  = 

1, 

IN+JN 

RL(I) 

:  x 

co-ordinate 

of 

trajectory 

at 

time 

T(I), 

I  = 

1,  IN+JN 

RM(I) 

:  y 

co-ordinate 

of 

traj  ectory 

at 

time 

T(I), 

I  = 

1,  IN+JN 

RN(I) 

:  z 

co-ordinate 

of 

trajectory 

at 

time 

T(I), 

I  = 

1,  IN+JN 

RMISS  :  miss  distance  parameters,  as  follows 


BMISS(l), 
RMISS (2) , 
RMISS (3) , 
RMISS (4), 
RMISS (5), 
RMISS (6), 
RMISS (7) , 
RMISS (8), 


time  of  nearest  miss 

x  co-ordinate  of  trajectory  at  time  of  nearest  miss 
y  co-ordinate  of  trajectory  at  time  of  nearest  miss 
z  co-ordinate  of  trajectory  at  time  of  nearest  miss 
range  at  time  of  nearest  miss 

1  component  of  trajectory  direction  cosine  vector 
m  component  of  trajectory  direction  cosine  vector 
n  component  of  trajectory  direction  cosine  vector 


//  JOB 

LOG  DRIVE 

CART  SPEC 

CART  AVAIL 

PHY  DRIVE 

0000 

0006 

0006 

0000 

V2  MO 9  ACTUAL  16K  CONFIG  16K. 

//  FOR 

♦  LIST  ALL 

♦  ONE  WORD  INTEGERS 
♦EXTENDED  PRECISION 
♦NAME  SINCA 

♦  ♦  SINGLE  CAMERA  SOLUTION  SUBROUTINE 
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SINGLE  CAMERA  SOLUTION  SUBROUTINE 
SUBROUT  I NE  S INCA ( T *RL*RM» RN » I N  » JN  »X • V * RMI SS • D ) 

DIMENSION  T  (  25  5  »RL ( 25  )  *RM (  25 )  *RN( 25 )  »X  (  3 )  tRNilSS  (  8  )  *D ( 3  *3 )  »A  l  9 )  » 
IS ( 25 ) 

C  CALCULATE  INITIAL  VALUES  FOR  (L*M*N) ♦  ROTATING  AXES  IF  NECESSARY 

RLL=RM ( 1 ) *RN ( IN )“RN ( 1 ) *RM (IN) 

RMM=RN < 1 ) *RL (IN) -RL ( 1 ) *RN (IN) 

RNN  =  RL ( 1 ) *RM  < IN ) -RM (1 ) *RL C IN ) 

RR=SQRT ( RLL*RLL+RMM*RMM+RNN*RNN ) 

RLL=RLL/RR 
RMM*RMM/RR 
RNN=RNN/RR 
K2  =  l 
J  =  1 

IF ( ABS(RNN)-O.l) 100*100*107 

100  IF ( ABS ( RMM )-0 *1)1 04 *104 *101 

101  KK  =  2 

102  DC  103  I  *  1 » 1 N 
RR=RN< I ) 

RN ( I ) =RM ( I  )*J 

103  RM( I )=-RR*J 
RR=RNN 
RNN=RMM*J 
RMM=-RR* J 

GO  TO  ( 108*109) *  K2 

104  KK  =  3 

105  DO  106  I  =  1  ♦  I N 
RR=RN( I ) 

RM ( I ) =RL ( I )*J 

106  RL ( I ) =-RR*J 
RR=RNN 
RNN=RLL*J 
RLL=-RR*J 

GO  TO  { 108*109) *K2 

107  KK  =  1 

108  I F ( RNN >10*11*11 

10  RLL=-RLL 

RMM=-RMM 

RNN=-RNN 

C  CALCULATE  COEFFICIENTS  FOR  THE  FIRST  ITERATION*TO  FIND  (L*M) 

11  DO  12  1=1*6 

12  A ( I ) =0*0 

DO  13  I  =  1 » I N 
A ( 1 ) =A ( 1 ) +RL ( I )*RL (  I ) 

A ( 2 ) =A ( 2 ) +RM ( I ) *RM (  I  ) 

A( 3 ) =A{ 3 ) +RN ( I )*RN (  I  ) 

A ( 4 ) =A ( 4 ) +RL ( I )*RM( I ) 

A ( 5 ) =A ( 5 ) +RN ( I >*RL{  I  ) 

13  A ( 6 ) =A { 6 ) +RM ( I >*RN( I  ) 

A ( 1 ) = A ( 1 ) -A { 3 ) 

A( 2 ) * A ( 2 ) “A ( 3 ) 

C  PERFORM  FIRST  I TERAT I  ON . AND  IMPROVE  ESTIMATE  FOR  (L»M) 

DO  14  1=1*30 
RLL2=RLL*RLL 

RMM2  =  RMM*RMM 
RNN2=RNN*RNN 
RR= 1 • 0/RNN 

F 1 =RLL*A ( 1 ) +RMM*A ( 4 ) +RR* ( l RNN2-RLL2 ) *A ( 5 ) -RMM*RLL*A { 6 ) ) 

F2  =  RMM#A  (  2  )  +RLL*A  (  4 )  +RR*  {  ('RNN2-RMM2  ) *A ( 6 )  -RMM*RLL*A  (  5  )  ) 
RR=RR/RNN2 

F3  =  A ( 4 ) +RR* ( ( RMM2-1 ) #RMM*A C  5 )  +  (  RLL2-1 ) *RLL#A ( 6 ) ) 

F4=A ( 1 ) +RR* ( (-3*RNN2-RLL2 )*RLL*A(5 ) + ( RMM2-1 ) *RMM*A ( 6 > ) 
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SINGLE  CAMERA  SOLUTION  SUBROUTINE 


F5  =  A(  2  1  +RR*  (  {  -3*RNN2-RMM2  )  *RMM*A  {  6  )  +  (  RLL2-1  >  *RLL*A  {  5  )  ) 
RR=F4*F5-F3*F3 
DM* ( F4*F2-F3*F1 ) /RR 
DL*(F1*F5-F3*F2)/RR 
'  RMK=RMM-DM 
RLL=RLL-DL 

RNN*SQRT( 1.0-RMM*RMM-RLL*RLL) 

RR  =  SGRT ( DM*DM+DL*DL  > 

IF(RR-O.COOOOOl) 15. 15*14 

14  CONTINUE 

15  J  — 1 
K2  =  2 

C  ROTATE  BACK;  TO  ORIGINAL  AXES  SYSTEM  IF  NECESSARY 

GO  TO  { 109  »  102  *  105  )  ♦« 

C  CALCULATE  MODIFIED  DIRECTION  COSINE  VECTORS  AND  THE  VECTOR 

109  DO  16  I  =  1 » IN 

RR=RLL*RL( I )+RMM*RM I ) +RNN*RN ( I ) 

KL  < I )*RL( I )-RR*RLL 
RM(  I  J  =  R M {  I  ) “RRfcRiMM 
RN( I )=RN( I )-RR*RNN 

RR=SORT(RL( I )*RL( I )+RM( I )*RM I )+RN( I )*RN< I ) ) 

RL C  I  )  *RL(  I  ) /RR 
R  M ( I ) ®RM ( I ) /RR 

16  R N ( I ) =  RN  ( I ) /RR 

RL ( 2  5 ) =  RM 11) *RNN-RN ( 1 ) *R.MM 
RM( 25 ) =  RN ( 1 ) *RLL”RL( 1 ) #RNN 
RN ( 2  5 ) «RL 1 1 ) *RMM-RM ( 1 ) *RLL 
C  CALCULATE  THE  VARIABLES  SHI 

DO  17  1  =  2. IN 

RR=RL( 1>#RLC I )+RM( 1)*RM( I ) +RN ( 1 ) *RN ( I > 

F4  =  RL 125) *RL ( I )+RM(25  >*RMC I  H-RN ( 25 ) *RN ( I ) 

17  S( I ) =RR/F4 

i*<  in-d/2 

DM= 1 • 0 / ( T ( I N ) -T (  1 )  ) 

DL= 1 • 0 / ( T ( I  +  l)-T( 1)  1 

ANN*  (  DL-DM  )  /  (  S  (  I  +  D-SUN)  ) 

AMM=A\N*S( IN) -DM 

C  PERFORM  SECOND  ITERATION  AND  IMPROVE  ESTIMATE  FOR  ( M » N ) 

DC  20  J  =  1 » 30 
DO  13  1=1.9 

18  AH  )=0.0 

DO  19  1=2. IN 

DL= 1 • 0/ ( ANN*S ( I J-AMM) 

DM=T ( I ) -T ( 1 ) -DL 
■  RMM2=DM*DM 
RLL2=DL*DL 
A ( 1 ) = A ( 1 ) +RMM2 
A  <  2 ) =A ( 2 ) -DM#RLL2 
A ( 3 ) = A ( 3 ) +S ( I >*RLL2*DM 
A{4) =A(4)+RLL2*RLL2 
A { 5 ) =A( 5)-2*DM*RLL2*DL 
A ( 6 ) = A ( 6 ) +S { I )*RLL2*RLL2 
A ( 7 ) =A ( 7 ) +2*S ( I )*DM*RLL2*DL 
A(8)=A(8)-2*S( I )  *S  (  I )*DL*DM*RLL2 

19  A(9)=AI9)+S1I)*S1I )*RLL2*RLL2 
RLL=AMM*AMM+ANN*ANN 
F1=AMM*A { 1 ) +RLL*A ( 2 ) 

F2=ANN*A ( 1 ) +RLL*A ( 3 ) 

F3=A { 1 ) +4*AMM#A ( 2 ) +RLL#A ( 4 ) +RLL*A { 5 ) 

F4=2* ( AMM*A ( 3 ) +ANN*A ( 2 ) ) +RLL# ( A < 7 ) -A ( 6 ) ) 
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F5*A( 1 ) +4* ANN* A I 3)+RLL*(A(8)+AI9) 1 
RR»1#0/(F4#F4“F5*F3) 

DKMF4*F2-F5*F1)*RR 

DL-(F4*F1-F3*F2)*RR 

AMM-AMM-DM 

ann»ann-dl 

RR-SORT  <  DM#DM+DL*DL ) 

IF ( RR-0«0000001 ) 21 *21 *20 

20  CONTINUE 

C  CALCULATE  TRAJECTORY  DATA  AND  MISS  DISTANCE  PARAMETERS 

21  DL*5QRT ( AMM*AMMtANN*ANN ) 

RR-V/DL 

OL-AMM#(TUN)-Ta)  )+X*0 
DM-ANNMTUN1-TI1)  ) 

RL I  IN) *RLI 25 )*DM+RLI 1 1 *DL 
RM I IN ) «RM I 2S )*DK+RM( 1 1  *DL 
RN  ( IN ) «RN ( 29  )*DM4>RN ( 1  )*DL 
RLll>-RR*RL(l)-Xm 
RM ( 1 ) *RR*RM ( 1 )"X( 2 ) 

RN ( 1 ) *RR*RN( 1 )»X( 3 ) 

RL( IN)«RR*RL( IN)-X(l) 

RM( IN)»RR*RM( INJ-XI2 ) 

RNI IN)«RR*RN( 1N)-X(3) 

DO  800  I  *  1  •  3 

A ( I ) ■RL( 1 )*D ( I «1 )+RM( L)*D{ I i2 l+RNI 1 )*D( I *3 ) 

800  A l  I  +  3)-RL(IN)*Dt  l.l)+RMUN)*Dl  l  •2H-RNUN)*DU  »3) 

F1*1.0/IT11N)-T<1) I 
RLL2*  t AI41-AI 1)  )*F.l 
RMM2* I  A I  5 ) -A 1 2 1 ) *F l 
RNN2* ( AI61-A (3 1 >  *F 1 
J*  I N+JN 
DO  23  l-l.J 
RR*T ( I 1-T 1 1 1 
RL(  I)  *A(  1 )+RR*RLL2 
RMU)»AI2)+RR*RMM2 
25  RN(I)-AC3H-RR*RNN2 

RR-- ( RL ( 1 ) *RLL2+RM ( 1 ) *RMM2*RN 1 1 )*RNN2 » /  C  V*V ) 

RMISS< 1)»RR+T(1) 

RMISSI 2)*RLI  1)*RR*RLL2 
RMISSI 3 ) «RMI 1 1 +RR*RMM2 
RM1SS ( 4) "RN 1 1 ) +RR*RNN2 

RMISS ( 5 ) ■SORT I RMISSI 2 1 *RMI$Sl 2)*RMI$S I 3)*RMISS( 3 )+RMISS(4 )*RMISS(  4 
1)) 

RM I SS 1 6 1 "RLL2/V 
RMISS  (  7 )  *RMM2/V 
RMISSI 8)*RNN2/V 
RETURN 
END 

VARIABLE  ALLOCATIONS 


AIR 

)» 0018-0000 

SIR 

)«QQ63-Q01B 

RLLIR 

1*0066 

RMM  ( R 

)  *0069 

RNMR 

1-006C 

RLL2IR 

)  *0072 

RMM2  ( R 

)  *0075 

RNN2IR 

1*0078 

F1IR 

1  *007B 

F2IR 

1-007E 

FAIR 

)  *0084 

F5IR 

1*0087 

DMIR 

1 *008A 

BUR 

)  -008D 

Ann  (  r 

1-0090 

<2  I  I 

)  *00A2 

J(I 

»*00A3 

KK 1 1 

) *00A4 

1(1 

1-00A5 

RRIR  )»QQ6F 
F3IR  1-QQBl 
AMMIR  ) *0093 


STATEMENT  ALLOCATIONS 


100 

■  0203 

101 

*0206 

102 

-Q2QF 

103 

11 

-02BF 

12 

•  02C3 

13 

•0318 

14 

20 

-07E3 

21 

•  07EC 

60U 

■  Q8E5 

25 

•0225 

104 

•0259 

105 

■02SD 

106 

•  0273 

107 

•02A7 

100 

•  02AB 

•0494 

•0977 

15 

■049D 

109 

•04AD 

16 

•053F 

17 

•05D3 

18 

■062F 

-02B0 

•OoDB 
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ONE  WORD  INTEGERS 
EXTENDED  PRECISION 


CALLED  SUBPROGRAMS 

ESQRT  EABS  EADD  EADDX  ESUB  ESUBX  EMPY  EMPYX  tDlV  ELO  ELOX  ESTO  ESTOX  ES0R 

EDVR  FLOAT  3U0SC  SNR  SUBIN 


REAL  CONSTANTS 

•100000000E  00*0004  •OOOOOOOOOE  00-0007  ilOOOOOOOOE  O1-QO0A  •1000Q0000E-06-00B0 
INTEGER  CONSTANTS 

l-OOCO  2-00C1  3-00C2  6-OOC3  30-0QC4  9«OOC5  4-OOC6 

CORE  REQUIREMENTS  FOR  SINCA 
COMMON  0  VARIABLES  180  PROGRAM  2408 

RELATIVE  ENTRY  POINT  ADDRESS  IS  00C7  IHEX) 


A  ppondlx  B 


ESBRX 


END  OF  COMPILATION 
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Fig. I 


Camera 


Fig. I  Camera  -  trajectory  plane 
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